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We explore the possibility that the observed onset of the Pioneer anomaly after Saturn 
encounter by Pioneer 11 is not necessarily due to mismodeling of solar radiation pressure 
but instead reflects a physically relevant characteristic of the anomaly itself. We employ 
the principles of a recently proposed cosmological model termed “the theory of inertial 
centers” along with an understanding of the fundamental assumptions taken by the Deep 
Space Network (DSN) to attempt to model this sudden onset. Due to an ambiguity that 
arises from the difference in the DSN definition of expected light-time with light-time 
according to the theory of inertial centers, we are forced to adopt a seemingly arbitrary 
convention to relate DSN-assumed clock-rates to physical clock-rates for this model. We 
offer a possible reason for adopting the convention employed in our analysis; however, we 
remain skeptical. Nevertheless, with this convention, one finds that this theory is able to 
replicate the previously reported Hubble-like behavior of the “clock acceleration” for the 
Pioneer anomaly as well as the sudden onset of the anomalous acceleration after Pioneer 
11 Saturn encounter. While oscillatory behavior with a yearly period is also predicted 
for the anomalous clock accelerations of both Pioneer 10 and Pioneer 11, the predicted 
amplitude is an order of magnitude too small when compared with that reported for 
Pioneer 10. 
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1. Introduction 

First reported in Ref. 1 and detailed in Ref. 2, the Pioneer anomaly refers to an 
unmodeled, blueshifted frequency drift with near uniform rate of change observed in 
the DSN residuals for both Pioneer 10 and Pioneer 11 spacecrafts with behavior in 
the latter portion of their cruise phases out of the Solar System resembling that of 
an additional near-constant acceleration in a yet to be unambiguously determined 
direction back toward the Solar System.^ While there have been many suggestions 
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and avenues of exploration for what could possibly account for this discrepancy,^’^ 
after a recent analysis of an extended set of DSN frequency data, there appears to 
be significant support for the origin of this anomaly as a recoil force associated with 
an anisotropic emission of thermal radiation off the spacecrafts due to the waste 
heat emitted by the radioisotope thermoelectric generators (RTGs) onboard. 
Nevertheless, the reported onset of the anomaly after Saturn encounter by Pioneer 
11 (see e.g. Refs. 13-15 and 16) remains an open source of concern after the same 
analysis of the extended data set appeared to be consistent with the possibility of 
this behavior, although taken to likely be the result of mismodeling of solar radia¬ 
tion pressure.Furthermore, other works independent of articles concerned with a 
thermal explanation effectively showed that the Pioneer anomaly as an extra physi¬ 
cal acceleration directed towards the Sun acting upon all Solar System bodies would 
be inconsistent with the observed motion of planets, spacecrafts, and other celestial 
bodies.^®’It would appear from these analyses that the Pioneer anomaly, 
if caused by physics not modeled by the DSN and not simply an unaccounted-for 
acceleration due to the intrinsic structure of the spacecrafts, would have to affect 
only the Doppler residuals calculated by the DSN and not the actual paths of our 
Solar System objects relative to one another. We will emphasize this point later 
in our discussion and argue that, under our proposed mechanism, objects in the 
Solar System will continue to orbit the Sun in the same manner as predicted by 
general relativity and expected by the DSN. While there have been several other 
documented anomalies in the Solar System that should be mentioned (reviewed in 
Refs. 33 and 34), such as the flyby anomaly,^® the secular increase in the astronom¬ 
ical unit (AU),^® an anomalous secular eccentricity variation in the Moon’s orbit, 
an alleged orbital anomaly reported in Saturn’s perhelion which seemed to disap¬ 
pear upon further analysis with Cassini data,^®“'*° and possible anomalous behavior 
reported in Jupiter’s perhelion albeit with a low level of statistical significance,^^’ 
our concern for the present article remains solely on the Pioneer anomaly. 

Our approach in this work is to explore the possibility that the reported on¬ 
set is a physical characteristic of the anomaly itself and thereby attempt to model 
this behavior near Saturn encounter along with the near-constant nature of the 
anomalous ‘acceleration’ long after for Pioneer 11. Furthermore, we look to model 
the near-constant nature of the anomalous ‘acceleration’ throughout the entire data 
set analysed for Pioneer 10. We assume throughout our discussion that, in future 
work when fitting to the extended set of DSN Doppler data, the thermal con¬ 
tribution to the anomaly can be approximately modeled according to a weighted 
mean of the original analysis of Ref. 2 with the analysis of Ref. 4, giving almost 
^ 12% of the full observed unmodeled ‘acceleration’ value, reported in Ref. 2 to 
be Up = (8.74 ± 1.33) x 10“^^ km s“^. We will discuss this assumption about the 
thermal contribution to the anomaly in more detail in our conclusions. To poten¬ 
tially accomplish this task of modeling the remaining nearly 8.0 x 10“^^ km s“^ 
of the anomalous ‘acceleration,’ we employ a recently proposed cosmological model 
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termed “the theory of inertial centers”and use the analysis of Refs. 44 and 45 
as inspiration for our reasoning. However, we must emphasize that we are not em¬ 
ploying the results of Refs. 44 and 45. While there have been other examinations 
of the possibility that the Pioneer anomaly can be linked to a cosmological origin 
such as the analyses in Refs. 46,47 and 48, we will not employ any of these results 
either given the assumptions of this newly proposed theory. Instead, with an un¬ 
derstanding of the fundamental assumptions taken by the DSN, we derive in this 
article our own results relevant to the predicted DSN frequency residuals accord¬ 
ing to the theory of inertial centers. Taking input data from the HORIZONS system 
of Jet Propulsion Laboratory (JPL) solely for calculating approximate values of 
our predicted residual term (i.e. our analysis is a rough ‘first-order’ approximation 
prior to proper fitting to DSN frequency data) and adopting a seemingly arbitrary 
convention due to a mathematical ambiguity to relate physical clock-rates in this 
model to DSN-assumed clock-rates, we discover that this theory is able to repli¬ 
cate the near-constant and blueshifted nature of the anomalous ‘acceleration’ for 
the latter interval of cruise phase out of the Solar System for both Pioneer 10 and 
Pioneer 11 as well as a sudden onset in anomalous ‘acceleration’ after Saturn en¬ 
counter by Pioneer 11. Thus, we tentatively suggest that the predictions of this 
cosmological model could provide an alternative to the thermal explanation of the 
Pioneer anomaly. However, a truly convincing explanation for the adoption of our 
particular convention for relating clock-rates is paramount to our understanding. 
Unfortunately, this resolution has yet to be determined. Still, given the ability with 
this convention to predict an onset for Pioneer 11 after Saturn encounter, we sug¬ 
gest the continuation of related work to further explore the consequences of these 
predictions in the form of potential observables for ongoing experiments in our solar 
system. 

This article is organized in the following manner. We begin with a brief review 
of the concepts employed in Ref. 44 and thoroughly detailed in Ref. 49. These are 
vital to our understanding of how the DSN calculates frequency residuals. We then 
review the principles presented in Ref. 43 most relevant for our purposes, subse¬ 
quently deriving an additional clock-drift term according to the theory of inertial 
centers for frequency residuals calculated by the DSN using the geodesic solutions 
of Appendix A. As mentioned previously, we employ a particular convention for the 
relation between physical clock-rates in this theory and DSN-assumed clock-rates in 
our derivation. We compare and contrast the “clock acceleration” (see section V.B 
of Ref. 2) predicted here with that of Ref. 44, emphasizing an additional weighting 
factor in the theory of inertial centers’ prediction dependent upon the direction of 
propagation of our transmitted photons relative to the direction of the center of 
our galaxy. This additional weighting factor, unique to the arguments presented in 
this article, produces a sudden onset of the predicted anomalous behavior after Sat¬ 
urn encounter by Pioneer 11. We provide concluding remarks and mention relevant 
future work. 
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2. Results and Discussion 

As alluded to in our introduction, the analysis in this section is very much inspired 
by the work of Ref. 44, which attempts to provide an explanation for the Pioneer 
anomaly within the context of our current standard model of cosmology. We offer a 
brief review to prepare for our own approach with the theory of inertial centers. The 
key feature pointed out in Ref. 44 is that current astrometric practice ignores the 
difference between coordinates suitable for an expanding universe and coordinates 
appropriate for globally flat space-time in general relativity after “factoring out” 
effects from gravitational sources. Presently, we take 

c(t2 - ti) = ^-21, (1) 

where r 2 i = |r 2 — ri| is the Euclidean coordinate distance between the point of 
emission and the point of reception of light. The DSN assumes this is equal to 
the proper distance ^21 in our empty vacuum scenario. Furthermore, c is the DSN 
constant for the speed of light in vacuum and f 2 ~ ti is the change in coordinate 
time between emission and reception. 

However, in the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric,one 
must take into account the rate of assumed expansion of the universe with the 
Hubble constant to linear order for our current experimental precision with 

photons. To account for this expansion. Ref. 44 suggests that the above equation 
must be replaced by 


c(A 2 — Ai) — r2i, (2) 

where 

= ti + —(3) 

analogous to the phenomenological “quadratic in time model” of Section XI.E in 
Ref. 2. The Newtonian approximation is assumed so that each observer is taken 
to be nearly stationary relative to an observer abiding by barycentric coordinate 
time (TCB) such that t = t, where r is the proper time of the observer measuring 
frequency shifts. In other words, it is suggested that to continue using a light-time 
relation of the form in (1) for empty space while still accounting for an expanding 
universe when modeling away from gravitational sources, one must take the light- 
ray coordinates to 0{Hq) to be {A,r} and not {t, r}. Ref. 44 goes on to find that 
the two-way fractional frequency shift calculated by the DSN for one round-trip will 
have an additional term 

^ = I-2[/521-J?o(t2-tl)], (4) 

1^1 

which would give the desired blueshift to possibly explain the anomalous acceler¬ 
ation in frequency shift seen with the Pioneer spacecrafts once we integrate over 
multiple round-trip cycles as the clock acceleration at quoted in Ref. 2 using the 
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formula 


2at{t2 — ti) 


( 5 ) 


observed model 


is similar in value to that of the Hubble constant. /32i is the relative velocity between 
the spacecraft and observer on Earth along the direction of motion of the light signal 


given by k 2 i • (v 2 — vi)/c and k 2 i is the unit vector directed toward the spacecraft 


from the observer on Earth. 

With the above summary as a preface, we apply this understanding of current 
practice to the predictions of the theory of inertial centers. This theory suggests that 
in order to account for the observed effects which our current standard cosmological 
model (i.e. A-CDM) attributes to the theoretical concepts termed ‘dark energy’ and 
‘dark matter,’ one need not assume the existence of these concepts, but instead one 
may reproduce this behavior by redefining what we consider to be inertial motion as 
well as what constitutes an inertial frame of reference in globally fiat space-time. 
However, locally within these redefined inertial systems, this theory was shown to 
abide by the axioms of general relativity.The conclusions presented in Ref. 43 
and differences with respect to general relativity result from the form of the affine 
parameter assumed in this theory for massive geodesics. 



( 6 ) 


where '/A is taken to be a universal constant and assumed to be our measured value 
of the Hubble constant. Furthermore, r^r) is the observer’s physical distance to the 


center of one of these newly defined inertial systems and r represents the physical 
clock-rate of the observer in this particular inertial system. Meaning, the physically 
observable elapsed time as measured by a clock carried along a given curve, denoted 
as proper time r, is not the affine parameter for the theory of inertial centers and 
thus is not invariant in this theory. Therefore, observable clock-rates depend upon 
both spatial position in a particular inertial frame as well as in which inertial frame 
the observer is observing. In contrast, the affine parameter for massive geodesics 
in general relativity takes the form X = c J dr, where instead c is taken to be a 
universal constant as the speed of light in globally flat space-time (see ch. 4.2 and 
4.3 of Ref. 53). 

To be explicit with regard to inertial motion, according to the theory of inertial 
centers, a massive object follows a geodesic path affinely parametrized by x (and 
not by r as in general relativity) when subjected to no net external forces with 
equations of motion given by 


0 = 


( 7 ) 


where 17“ is the tangent vector to the geodesic curve denoted x° (x) in a particular 
coordinate system (i.e. [/“ = dx°‘/dx). Notice that we have employed abstract index 
notation and will continue to do so throughout this article (see ch. 2.4 of Ref. 53). 
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Expanding (7) in a particular coordinate system, one finds that the equations of 
motion for our massive object reduce to 


dx^ 


where r“j^ are the components of the affine connection for the coordinate chart that 
we are employing to describe the geodesic path. However, to describe this path in 
terms of physical observables (i.e. proper time r), one must use the conjectured 
relation between x and proper time in the particular inertial system within which 
the massive object moves. We have 


0 = 


dx°‘ 

fdrV 

d^a:“ dx^ dx'' 

dx^ dr 

\dx) 

dr'^ ' dr dr 


(9) 


where dx/dr = \/A • r(T). Additionally, if the bracketed term in (9) were to equal 
zero, it is evident that we would have the geodesic equations of motion according 
to general relativity,®^ and thus (9) shows the clear physical difference between 
the geodesics of the theory of inertial centers and those of general relativity (i.e. 
in general, d^r/dx^ ■ dx°^/dT ^ 0; the general relativistic limit occurs for ^(t) = 
constant). In terms of the proper velocity of the massive object with components 
= dx°‘/dr, our geodesic equation implies 

1 fl'P 

= -—u\ ( 10 ) 

r dr 

in accordance with eqn. (3.3.2) of Ref. 53. 

To describe one of these newly proposed inertial frames of reference, it appears 
“most natural” to employ the radial Rindler chart 

ds^ = —Ar^dt^ + dr^ + cosh^{'/At){d9^ + sin^ ddcj)^), (11) 

where = Hq is taken to be the Hubble constant and = —dx^ = —Ar^dr^ 
for massive geodesics, as coordinate time t in this chart corresponds to the physical 
time elapsed on the clock of a stationary observer according to the assumptions of 
this theory (i.e. inertial time). Notice that observable clock-rates in this theory differ 
from those of general relativity due to the form of the affine parameter and reflected 
in the fact that, in globally flat space-time, they are dependent upon both the relative 
velocity between observers as well as relative spatial position with respect to the 
center of the inertial frame. Ref. 43 terms the center of each new inertial frame an 
“inertial center” (in this chart, r = 0 is taken to be the location of our inertial 
center). 

In the following, we use the term “stationary observer” to refer to any observer 
that is stationary relative to the center of the inertial frame in question as defined 
according to the theory of inertial centers. Furthermore, Ref. 43 takes the Milky Way 
itself to constitute one of these new inertial frames with the center of the Milky Way 
in the direction of Sagittarius A* constituting the center of the inertial reference 
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frame relevant for observers in our solar system. In other words, it is assumed in this 
model that at the center of our galaxy there exists an inertial center and not a super- 
massive black hole (for recent tension with our current understanding of black hole 
formation, see Refs. 54 and 55). However, as Ref. 44 mentions in his work, the adop¬ 
tion for the time coordinate in current programs assumes the Minkowski coordinate 
parametrization of space-time to be accurate for completely empty space, yet this is 
also not the case for the theory of inertial centers. Meaning, the DSN assumes that 
Minkowski time is consistent with the progression of a physical clock in deep-space 
far from gravitational sources, which in fact is not the case when operating under 
the assumptions of this theory. Therefore, we should formulate calculated frequency 
shifts in terms of the coordinate time assumed by the DSN to determine possible 
residuals. We emphasize the word “calculating” because this additional ‘frequency 
shift’ is not actually a physical shift in the frequency of our two-way electromagnetic 
signals according to this model (see the end of Appendix B). Instead, it is conjec¬ 
tured that the DSN’s definition of time resulting from the assumed validity of (1) 
in empty space is inconsistent with the manner in which physical clocks progress. 
This approach is in contrast to the naive procedure outlined in “Application to a 
local gravitational system” in Ref. 43 where the Pioneer anomaly is suggested to be 
the result of a physical frequency shift. 


2.1. Stationary observer treatment 


For the general case of examining the two-way frequency residuals as calculated 
by the DSN from light signals sent between a spacecraft in our solar system and 
our tracking stations here on Earth, we must analyse the range-rate calculation 
employed in Ref. 49 with the understanding that the DSN’s definition of time does 
not coincide with physical clock-time according to the theory of inertial centers. 
Thus, we look to relate clock-rates by deriving the differential change in DSN- 
assumed coordinate time with respect to physical time in this theory. Our approach 
to accomplishing this task relies upon relating the light-time expected by the DSN 
with the light-time expected by the theory of inertial centers as this quantity, in 
essence, defines changes in time for a fixed proper distance traversed (i.e. a light 
clock). We refer to the geodesic solutions for a light signal traveling within the Milky 
Way inertial system with an inertial center assumed to be located in the direction 
of Sagittarius A*. From (A.22) and (A.26) in Appendix A, we have for the radial 
solution of the geodesic path of a massless particle 


r{t) = ro 


^ sinh(VAt) -t 
Co 


1 + 



2 

sinh^(-\/At) 


( 12 ) 


where Vq = dr/dt\t=o is the initial radial velocity of our particle at < = 0, tq = 
r{t = 0) is the initial distance to our inertial center, and cq = Additionally, 

Vo = Co = y/'^i v'^Vi\t=o is the initial speed of our massless particle where this 
sum occurs over a constant t spatial hypersurface of our metric (11) with induced 
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Riemannian metric components hij. Solving for the elapsed coordinate time t for a 
light signal traveling between rg and r 


t = —= sinh ^ 

Va 


^ sinh ( ln(r/r-o) 


L^o 


(13) 


but rij = co(ko-ro): where ko is the initial direction of propagation of our light 
signal and fg is the initial radial unit vector for the emitter of said signal within 
our inertial system. Then, 


t± = —sinh ^ 

Va 

where is the elapsed coordinate time for a signal traveling away from our inertial 
center and is the elapsed coordinate time for a signal traveling toward our inertial 
center. For a two-way signal sent between nearly stationary observers labeled 1 and 
2 within the galactic inertial system, notice that 


Ikn-fnl 


sinh ±ln(r/ro) 


(14) 


^uplink — sinh 
^downlink — 7= sinh 

VA 


-1 


|krfi| 


sinh ±ln(r 2 /ri) 


-J—sinh ( Tln(ri/r 2 ) ) 

k2-r2 V / 


(15) 

(16) 


For signals sent within the Solar System, we will make the approximation ki jupUnk = 
k 2 1 uplink = —k 2 1 downlink and assume our galactic radial unit vector is approximately 
the same everywhere within our solar system (i.e. fi = f 2 = f). Thus, tupiink = 
idowniink- Then, the elapsed coordinate time for a light signal sent one-way between 
observers located at two radial distances rf, and Vg is given by 


t = —= sinh 

Va 


|k,-f| 


sinh ( \n{n/rg) 


(17) 


where (i.e. = r 2 , Vg = ri for r 2 > ri and ri, = ri, rg = r 2 for r 2 < ri). 

The DSN assumes the elapsed coordinate time Tdsn for light traveling a total 
physical distance d 2 i (i.e. proper distance) is given by^® 


Tdsn = —(18) 

CDSN 

where cdsn is the DSN-assumed value for the speed of light in empty vacuum. Yet, 
given our assumptions of the near-constancy of the direction of our galactic radial 
unit vector within the Solar System, we find that the DSN assumes the physical 
radial distance traversed by a one-way light signal in the galactic reference frame 
should be given in terms of DSN coordinate time by 


|Ar| Ri |ks-fI cdsnTdsn- 


(19) 


We will use this as our analogue to the defining relation for DSN light-time (18) in 
the rest of our analysis. As well, we take our proper radial length in the galactic 
frame according to the theory of inertial centers to be |Ar| = |r 2 — ri| = r^ — rg. 
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In looking for an expression relating DSN-assumed clock-rates to physical clock- 
rates in the theory of inertial centers (i.e. /dt) ^ we come upon the unfortunate 

consequence that our light-time in this model relies upon fractional differences in 
distance, whereas the DSN definition of light-time relies upon the actual difference in 
coordinate values. Thus, the question becomes: where do we place the dependence on 
Tdsn in our fractional expression rj, / Ts given that the DSN defines time based off of 
rij — rgl In other words, expressing the fractional difference in terms of | Arj = rt, —r^, 
one has 


^ _ 1 _l_ |Ar| _ rb 
Ts Ts rb-\Ar\' 


( 20 ) 


But the remaining presence of radial coordinate values rb or in this fractional 
difference expression, which we will eventually plug back into (17), inherently forces 
us to make a choice as to whether or not rb or has dependence upon Tdsn as we 
are looking for the differential change in range with respect to time elapsed on the 
physical clocks of this model. We focus on two possibilities: 


rb = rbiTusN) Ts = constant. 


( 21 ) 


or 


rb = constant = rs(rDSN), 

such that (19) reproduces a DSN-assumed null geodesic of the form 

?'b(TDSN) = |ks-f IcdSnTdSN + rs, 


( 22 ) 

(23) 


for the former case or 


rsiTusn) = — Iks-rjcDSNTbsN +rb, 


(24) 


for the latter. Again, we’ve assumed |ks-f| = |k(,-f| for the one-way signal. Note 
that these DSN-assumed null geodesic equations are valid only for our assumptions 
of an electromagnetic signal traveling within a small region of the galaxy. 

We seek to express our fractional difference in radial positions in terms of only 
assumed constants and Tdsn- For the former case, one has in (17) the expression 


In 



In 



= In 



|ks-r|cDSNTDSN 


(25) 


where we’ve used defining relation (19). If one were to adopt Tdsn as a time coor¬ 
dinate covering our entire inertial system, he/she would find that this expression is 
valid for all values of 0 < Tdsn < oo. However, for the latter case, one finds 


In 
= In 


= In 


rb 


rb - |Ar| 


rb 


rb — |ks-f IcdsnTdsn 


(26) 
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where this expression restricts our Tdsn ‘independent’ parameter values such that 
0 < TbsN < |cdsn- In reality, of course, Tdsn is not an independent param¬ 

eter as its measured value depends upon the proper distance between our two ob¬ 
servers. Nevertheless, referring to (19), one hnds, as of course we should expect, that 
this restriction on Tdsn is consistent with our expectations for the DSN-assumed 
light-time since we must have rg > 0 for our current analysis. Still, in this instance, 
it does not seem that one could use Tdsn as a global ‘time coordinate’ for a DSN 
parametrization of empty space-time in the inertial system given the restrictions 
imposed by the theory of inertial centers. While it seems rather arbitrary to suggest 
that this has any physical relevance on our interpretation and predictions for DSN 
observations as both choices for Tdsn are consistent with DSN-assumed light-time 
conditions, we mention the idea for discussion. Do the restrictions on the allowed 
values of Tdsn for the latter case point to the suggestion that the DSN time coor¬ 
dinate in this scenario is unphysical and thus the only physical solution lies with 
the convention rj, = ^^(Tdsn)? Given these arguments, we’ll initially approach this 
problem of determining DSN residuals according to the theory of inertial centers by 
adopting the former convention of rj, = rt,(TDSN), where Tdsn as an ‘independent’ 
parameter can take all values greater than zero in its relationship with t. 

Manipulating (17) to isolate our fractional difference factor, one finds 

^ _ gSinh“^[|ks.f| sinh(%/At)] (27) 

But expressing in terms of Tdsn, "we have for the former case the relation 


Tdsn 


Vs 

|kg-f |cdsn 


gSinh ^[|ks-fI sinh(\/Ai)] 



(28) 


where we again emphasize that the smaller radial coordinate value, r^, is taken to 
be a constant. Proceeding to the relevant observable, we find that the theory of 
inertial centers prediction for the differential change in DSN-expected light-time 
with respect to the physical clock of an observer stationary relative to the center of 
the Milky Way is given by 


dTDSN _ \/Arge™‘’ '[Ik.-fDinh(VAt)] cosh(yAt) 
~ "°SN ^l + (k^.f )2 3inh2(VAt) 


Expanding to first-order in vGVt, we find 


rfTpSN 

dt 


CDSN 


1 + |kg-f l-x/At 


+ 0{At'^). 


(30) 


However, the frequency shift calculated hy the DSN is taken to be the ratio of 
differential proper times, where it is assumed that the total light-time is given 
by Tdsn.^® Thus, we take the time coordinate employed by the DSN to be Tdsn 
and not t when calculating range-rate. Meaning, in this model, the two-way DSN 
fractional frequency shift [r'a/zzijDSN does not represent an accurate calculation of 
the fractional shift in the frequency observable itself (i.e. v^/vx 7 ^ [r'a/jzijDSN) due 
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to the fact that the DSN’s definition of time drifts away from actual physical time 
according to the theory of inertial centers with our first-order relation modeling this 
drift given by (30). 

Then, the two-way DSN fractional frequency shift between stationary observers 
{dt/drli = 1 and dTDSNli/dTbsNlj = 1 for i,j = 1,2,3) is predicted by the theory 
of inertial centers to be 




dTbsN 

f ^Tdsn 



dTbsN 

/dTbsN 


.^1. 

DSN 

dt 

2 V 

J J 

uplink 

dt 

3 V 

2 ) \ 


downlink 


1 -b |k2-f|VAt2 ^ 1-b |k3-r|VAf3 
1 -I- |ki-f|A/Atl 1-I- |k2-f|A/At2 
= 1 + |k-r|VA(f3 - h) + 


assuming ki = k2|upiink = — k2|downiink = — k3 = k. However, as we found earlier, 
tupiink = tdowniink for our assumptions in this small region of our galaxy, where 
tupiink — ^2 and tdowniink — ^3 ^ 2 - ^Ve have 


= l + 2\k-r\VA{t2-ti) + 0 {At'^), (32) 

DSN 

where t 2 — ti is the one-way light-time measured by our stationary physical clock. 
One finds for the two-way frequency shift as calculated by the DSN an additional 
drift term with clock acceleration given by 

at = -/Alk-fl. (33) 




Furthermore, with the adopted convention of ri, = rb(TDSN), our clock-drift term 
produces a blueshifted residual independent of whether or not the signal initially 
travels away from or toward the galactic center. Notice that one obtains the result 
of Ref. 44 for photons traveling either directly toward or away from the galactic 
center as a/A is taken to have the value of the Hubble constant in the theory of 
inertial centers. Since Pioneer 10 and Pioneer 11 traveled nearly directly away from 
and toward the center of our galaxy, respectively, on their way out of the Solar 
System^ (i.e. |k-f| k, 1 ), these results appear at first glance to be consistent with 
the late-time observations from Ref. 2 of a blueshifted clock-drift residual for both 
Pioneer 10 and Pioneer 11 as the observed at was similar in value to that of the 
Hubble constant. 

Nevertheless, the key to comparing these predictions with actual DSN residuals 
lies in calculating the additional drift term from (32). In order to proceed, we 
require the complete history of the path of the spacecraft relative to Earth’s path 
over the course of our observations to calculate the direction of propagation relative 
to the center of the Milky Way of each light signal sent between our stations on 
Earth and the spacecraft. Still, the above equation is very much an approximation 
because we assumed both the spacecraft and Earth observer are stationary relative 
to one another and chose to neglect gravitational effects. However, one can think 
of this analysis as “factoring out” the Doppler shift due to the relative velocity of 
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our observers as well as any potential frequency shift due to gravitational sources 
(see Appendix B for a more general derivation taking into account the Doppler 
effect). Meaning, additional terms predicted by this theory with respect to the DSN’s 
expectations in a full relativistic derivation will be suppressed by the multiplication 
of -x/At and thus can likely be disregarded to first-order in y/Xt, (v/c), and GM / c^dM 
where dM is the distance to the center of a massive object in our solar system. 
Furthermore, the weighting factor |k-f | in our clock acceleration reflects the amount 
that the DSN’s definition of time coincides with physical time according to the 
theory of inertial centers (see (30)), as the difference in light-time between that 
assumed by the DSN in the Solar System and what physically occurs in the theory 
of inertial centers is smallest when both observers are located the same distance 
away from an inertial center (i.e. |k-f| = 0) and largest when the photon travels 
radially between observers located at two different distances from the center of our 
galaxy (i.e. |k-f| = 1). 

While examining alternatives to the thermal explanation for the Pioneer 
anomaly, one may ask and rightly so: if the initial results in our DSN calcula¬ 
tions of two-way frequency shifts for Pioneers 10 and 11 are nearly the same as 
those produced by Ref. 44, how can we tell whether or not the Pioneer anomaly 
points to a consequence of the assumptions of the theory of inertial centers or to 
a consequence of the analysis of Ref. 44? Unfortunately, due to the paths of both 
Pioneer spacecrafts almost directly toward and away from Sagittarius A* on their 
way out of the Solar System, it would be difficult to tell the difference between the 
predictions of both models with the precision of the available data for the average 
late-time behavior of the anomaly. However, the results of Ref. 44 o,re based upon 
assumptions of isotropy. Thus, they should hold for spacecrafts moving in any par¬ 
ticular direction in outer space far from the Sun in the galactic frame. The same is 
not true for the theory of inertial centers as one must consider the paths of these 
spacecrafts relative to the center of the Milky Way. Therefore, if we were to send 
a spacecraft out moving tangentially to the direction of the center of the Milky 
Way, we would obtain a different result while working under the assumptions of 
the theory of inertial centers versus those of Ref. 44. For context in the rest of our 
analysis, the predicted residuals calculated by the DSN should vary according to 
the trajectory of the spacecraft relative to our DSN stations on Earth if one models 
with the theory of inertial centers. Thus, a gravitational assist that signihcantly 
alters the trajectory of a spacecraft, such as that of Pioneer 11 by Saturn, could 
potentially “turn on” the anomaly according to this theory, whereas this type of 
“turn on” would not seem to be possible according to the analysis of Ref. 44. 

Before proceeding to a more comprehensive analysis of our predictions for several 
spacecrafts, we must return to the latter scenario of our convention choices with 
fs = *'s(2Usn) for our dependence on DSN-defined time. We recall the restrictions 
on this ‘time coordinate’ 

0 < Tdsn < —yy-, 

|ks-r|cDSN 


(34) 
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which led us to the suggestion that we can possibly consider this scenario to be 
unphysical as this candidate for a DSN time coordinate cannot cover all relevant 
physical values in a relation between Tdsn and t. Still, as mentioned in our abstract, 
we are not completely convinced by this argument. Working through a similar anal¬ 
ysis as that employed for the former case, one finds to first-order in VM 


dt 


cdsn 


1 — Iks-flv^t 


OiAf), 


(35) 


which would ultimately produce a predicted redshift in DSN residuals. Therein lies 
the major issue. The adoption of a particular convention results in a change in sign 
of our drift term. Yet, we have chosen the former convention mainly because it is 
consistent with the observations from Ref. 2. Future work will require investigating 
the validity of our proposed argument that the unphysical nature of the restrictions 
imposed by the theory of inertial centers rules out the latter scenario. 

For emphasis in our next section when speaking of the consistency of these 
results with our current experimental precision of the ephemerides which tend 
to be dominated by ranging measurements,^®^®® we examine how the predictions 
of the theory of inertial centers affect DSN calculations of range. The important 
point to take away is that range represents a time-of-flight calculation for a sin¬ 
gle round-trip of light signals sent back-and-forth to objects in the Solar System. 
Contrast this with DSN Doppler observables (i.e. range-rate), which represent a 
change in Doppler cycle count over a predetermined amount of time where the 
cycle count is calculated from an integration of the cycles of frequency accumu¬ 
lated over the count time, thus producing a compounding effect when analyzed over 
long data intervals (see ch. 13.3 of Ref. 49). Theoretically, one can approximate 
the manner in which the DSN reproduces this compounding effect by multiplying 
together multiple fractional frequency ratios as if the signals are bouncing back- 
and-forth between Earth DSN stations and spacecrafts within the Solar System 
(i.e. Vnfvi = {vn!Vn-i){vn-i!Vn- 2 )''' (^ 2 1vi))■ However, this would technically be 
inaccurate when comparing with the actual implementation detailed in Ref. 49, as 
remarked upon in section IX.C of Ref. 44, because, after a signal’s two-way round- 
trip, the DSN retransmits another signal at the same frequency at which it initially 
emitted, vi, but time-stamps the point of reception of the previous signal (i.e. ts) 
and records the difference in frequency between the emitted and received frequency, 
1/3 versus vi. While one still accumulates phase in a “compounding” manner with 
the DSN-employed time-integration process, theoretically it appears simpler to see 
this effect with a multiplication of multiple fractional frequency ratios. Neverthe¬ 
less, our point here is to emphasize that ranging measurements do not have this 
“compounding” effect as they only measure a single round-trip for each light signal. 

To determine possible additional terms in DSN-range calculations according to 
the theory of inertial centers, we return to our expression for the DSN-calculated 
one-way light-time (28) and realize from (18) and (19) that given the manner in 
which our DSN time coordinate is defined, the actual physical range predicted by 
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the theory of inertial centers will he p = cdsn^dsn- Thus, we have for our range 
expression when employing our former convention (23) 


p = cdsnTdsn 


Ts 

|ks-r| 


gSinh '■[|ks-f I sinh(-\/At)] _ ^ 


(36) 


For light-times within our solar system, we examine terms to second-order in y = 
'/At by Taylor expanding about y = 0. One has 


which gives 


p{y = 0) = 0 
dp 
dy 


d?p 

dy'^ 


v=o 


= rs|ks-f|, 


y=o 


p = p{y = 0) + ^ 

dy 


y=o 

= y/Avs 


1 d^p 
2dy‘- 
1 , 


y^ + ... 


y=o 


t -I- -\ks-r\'/A/ + ... 


(37) 

(38) 

(39) 


(40) 


where t is the actual time-of-flight in this theory. As the DSN expects range to be 
given by"^® pdsn = cdsn • we’ll have as residuals ApDSN 


A/Odsn = Pdsn — P = (cdsn — VAr/t — y/Ar, 

' y/Kr, 


i|ks-f|A/At® -I-... 


= ( 1 - ^pdsn - i\/A|ks-f| 

Cdsn / ^ 


Cdsn 


f Pdsn \ 
\Cdsn / 


(41) 


But from “Reduction to special relativity” in Ref. 43 as well as Appendix A 
and Appendix B, one realizes that y/Arg actually represents the local speed of light 
that the observer located a distance from the center of our galaxy measures in 
this theory. Thus, the value adopted by the DSN for cdsn should correspond to a 
measurement for the local speed of light by a observer stationary relative to the 
galactic center and located within our solar system according to this model (in 
particular, measured at a point in Earth’s orbit about the Sun). Taking cdsn = 
■\/ArDSN, one can express Vg in terms of our DSN observer’s distance to the galactic 
center with 


fs = Cdsn + Ar, (42) 

where Ar is the difference in distance to our galactic center between the locations 
of our two hypothetical observers Vg and r^sN (he. range in the radial direction in 
the galactic frame). Notice that if we are thinking in terms of our DSN stations on 
Earth sending light signals out to our spacecrafts then r^/rDSN will only deviate 
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away from one for spacecrafts located closer to the galactic center and located at 
large distances away from Earth. One has 

ApDSN = -^^PDSN - lVA\k-r \(1 + +.... (43) 

However, for order of magnitude purposes, we examine the largest possible Ar value 
such that the spacecraft is in exactly the same direction as our galactic center from 
the Earth observer’s perspective. Then for Earth-Mars ranging and Lunar Laser 
Ranging (LLR) one expects as upper limits ApDSN -1 cm and A/Odsn ~ 10“^ cm, 
respectively, similar in magnitude to the range results obtained in Kopeikin’s anal¬ 
ysis (see section IX.B of Ref. 44). Nevertheless, our current experimental accuracy 
for ranging to Mars is of the order of ~ 10 m and for LLR ~ 1 cm. Therefore, these 
ranging experiments would be unable to detect the predictions of this theory. 


2.2. Predictions using data from HORIZONS 

To facilitate comparison with plots and values from Ref. 2, we formulate all of our 
predictions for this model in terms of an ‘acceleration’ on our spacecrafts, defining 
Op = CDSNOt = -x/AcDSNlk-fI, where cdsn is the DSN-assumed value for the speed 
of light in vacuum and, as mentioned earlier, we only investigate for the adopted 
convention of Vb = rb(TDSN) which predicts a blueshifted clock-drift residual. How¬ 
ever, one must keep in mind that this ‘acceleration’ is fictitious in the sense that the 
spacecraft is not experiencing an additional acceleration according to the theory of 
inertial centers. Instead this ‘acceleration’ is a reflection of the difference in clock- 
rates between physical clocks in the theory of inertial centers and those assumed 
by the DSN as the clock acceleration at is what is truly physically relevant for this 
model affecting only the DSN-calculated frequency data. Furthermore, in this the¬ 
ory, one would expect all objects in the Solar System to be moving along a similar 
geodesic path about the center of our galaxy and thus, while that geodesic should 
differ from what one would expect according to general relativity when analyzing in 
the galactic frame, relative accelerations between objects in the Solar System frame 
of reference should not differ significantly from the expectations of the DSN as this 
theory reduces to general relativity within a localized region of the galaxy when objects 
remain nearly the same distance away from the galactic center (see “Reduction to 
special relativity” and “Limitations of the study, open questions, and future work” 
in Ref. 43). Then, one would not expect the positions of objects in the Solar System 
relative to one another to differ significantly from DSN expectations when operating 
under the assumptions of the theory of inertial centers, and thus we do not expect 
astrometric observations in the ephemerides to be different from general relativis¬ 
tic expectations. As with Ref. 44, the only observable differences with our current 
experimental sensitivity should occur with light, and in particular Doppler mea¬ 
surements due to their “compounded” nature accentuating that difference, thereby 
leading to our potential explanation of the Pioneer anomaly. Nevertheless, as shown 
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at the end of our previous section, ranging measurements in the ephemerides will 
not be affected according to this theory given our current experimental precision. 

Since the majority of observations that comprise the ephemerides are ranging and 
astrometric/optical in nature, it does not appear that terms predicted by this theory 
would alter the current general relativistic expectations of the ephemerides. Thus, 
the predictions of this theory should not conflict with our current understanding 
of the motions of large massive bodies in the Solar System, experimentally verified 
to a very precise degree.^®’ 18-32,41,42,56-58 p^j. purposes, we approximate the 
DSN-assumed value for the speed of light in vacuum as cdsn = 3.00 x 10® km s“^. 

Table 1. Replication of table I from the original Pioneer anomaly analysis.^ Determinations of ap in units of 10“13 km 
from the three time intervals of Pioneer 10 data and from Pioneer 11. JPL analysis defined the time intervals as I (3 
January 1987 to 17 July 1990); II (17 July 1990 to 12 July 1992); and III (12 July 1992 to 22 July 1998). Results from 
various ODP-SIGMA and CHASMP calculations are listed. “WLS” signifies a weighted least-squares calculation. “BSF” 
signifies a batch-sequential filter calculation. “With solar corona model” designations refer to methods using the Cassini 
solar corona model. Lastly, “Corona, weighting, and F10.7” designation refers to using the Cassini solar corona model with 
corona data weighting and F10.7 time variation calibration. Errors given are only formal calculational errors. 


Program-Estimation method 

Pioneer 10 (I) 

Pioneer 10 (II) 

Pioneer 10 (HI) 

Pioneer 11 

SIGMA, WLS, no solar corona model 

8.02 

± 

0.01 

8.65 

± 

0.01 

7.83 

± 

0.01 

8.46 

± 

0.04 

SIGMA, WLS, with solar corona model 

8.00 

± 

0.01 

8.66 

± 

0.01 

7.84 

± 

0.01 

8.44 

± 

0.04 

SIGMA, BSF, 1-day batch, with solar corona model 

7.82 

± 

0.29 

8.16 

± 

0.40 

7.59 

± 

0.22 

8.49 

± 

0.33 

CHASMP, WLS, no solar corona model 

8.25 

± 

0.02 

8.86 

± 

0.02 

7.85 

± 

0.01 

8.71 

± 

0.03 

CHASMP, WLS, with solar corona model 

8.22 

± 

0.02 

8.89 

± 

0.02 

7.92 

± 

0.01 

8.69 

± 

0.03 

CHASMP, WLS, with corona, weighting, and F10.7 

8.25 

± 

0.03 

8.90 

± 

0.03 

7.91 

± 

0.01 

8.91 

± 

0.04 


We use a value near that quoted in table I of Ref. 2 for the anomalous ‘accelera¬ 
tion’ from interval III of the JPL analysis (12 July 1992 to 22 July 1998) of Pioneer 
10 data to determine an approximate value for i/A in this theory since the solar 
radiation pressure on the spacecraft is most diminished in this interval (i.e. Pioneer 
10 is farthest from the Sun in interval III; see fig. 6 of Ref. 2 for solar radiation 
pressure levels fitted during Earth-Jupiter cruise phase and Table 1 in this work for 
a replication of table I of Ref. 2). We take a/A = 2.67 x 10~^® s“^ for the Hubble 
constant in our current analysis corresponding to an interval III anomalous ‘accel¬ 
eration’ for Pioneer 10 of approximately Oplpio = 7.84 x 10“^® km s“^ (eqn. (23) 
of Ref. 2). We wish to emphasize that this article is concerned with the predicted 
behavior of the anomalous residuals after Saturn encounter for Pioneer 11 such that 
our focus is not primarily on the exact value of the anomalous ‘acceleration’ seen by 
both spacecrafts. Meaning, we have used input data generated from HORIZONS to 
determine an approximate direction of propagation of these electromagnetic signals 
(i.e. to determine our |k-f | term) solely for the purposes of a rough calculation of 
residuals predicted by the theory of inertial centers. Nevertheless, we keep in mind 
that the Pioneer anomaly provides an avenue for determining the value of the Hub¬ 
ble constant within our own galaxy according to this theory. Fitting to the extended 
set of DSN frequency data in future work will help us to appropriately revise this 
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value for '/A- 




Fig. 1. Side-by-side plots of the anomalous ‘acceleration’ predicted by the theory of inertial 
centers for Pioneer 11 assuming "s/A = 2.67 x 10“^® s“^ (a) and the trajectory of the spacecraft 
in terms of heliocentric distance (b) from 1 January 1977 to 1 January 1989. Note that Saturn 
encounter occurred on 1 September 1979,^^ subsequently resulting in a jump in the predicted 
anomalous ‘acceleration’ as can be seen in (a) near a heliocentric distance of 10 AU. Input data 
used to generate these plots was taken from the JPL HORIZONS system. 


Turning our attention to Fig. 1(a) and Fig. 1(b) where we have plotted the 
anomalous ‘acceleration’ Op for Pioneer 11 according to the theory of inertial centers 
versus heliocentric distance, we realize by comparing with fig. 7 of Ref. 2 that the 
initial results from the predictions of this theory are promising. One finds near a 
heliocentric distance of 10 AU after encounter with Saturn by Pioneer 11a jump 
in anomalous ‘acceleration’ ultimately plateauing into a near-constant Hubble-like 
value for the clock acceleration as was observed and reported. Furthermore, we also 
find oscillations in the anomalous ‘acceleration’ for the late-time behavior during 
plateau. We see from Fig. 2, where we’ve plotted the predicted Up for Pioneer 10 
versus date, that the theory of inertial centers prediction of a near-constant value for 
Op during the entire data interval plotted, where Pioneer 10 traverses a heliocentric 
distance range of approximately 20 to 70 AU, is consistent with the conclusions of 
section V.A in Ref. 2 where the observed anomalous ‘acceleration’ for Pioneer 10 
was found to be within a range of ±2 x 10“^^ km s“^ of (8.09±0.20) x 10“^^ km s“^ 
for a spacecraft heliocentric distance between 40 and 60 AU. We also find oscillatory 
behavior in Up for Pioneer 10 as well with the oscillations suppressed over time. 

It seems at first glance that this theory could provide a potential avenue to mod¬ 
eling the late-time annual oscillatory behavior discussed in section IX.C of Ref. 2 
as the near-annual period and decaying nature of the oscillations predicted for Pi¬ 
oneer 10 and Pioneer 11 in Fig. 1(a) and Fig. 2, respectively, are similar to that 
reported for Pioneer 10. Physically, one would provide the following reasoning be- 
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Fig. 2. Anomalous ‘acceleration’ predicted by the theory of inertial centers for Pioneer 10 from 
1 January 1981 to 22 July 1998. 


hind these oscillations. As the trajectories of these spacecrafts were relatively radial 
in the galactic system during their cruise phases out of the Solar System, we would 
interpret the late-time oscillatory behavior of the anomalous ‘acceleration’ seen by 
each as being due to Earth’s revolution about the Sun each year, since this revolu¬ 
tion produces an oscillatory change in the direction of propagation of our emitted 
photons relative to the galactic center (i.e. the factor |k-f | oscillates throughout the 
year due to Earth’s changing position in the Solar System). This, in turn, produces 
an oscillatory change in the amount the DSN’s definition of time overlaps with 
physical clock-rates according to the theory of inertial centers. Furthermore, this 
predicted oscillation is suppressed over time since k essentially becomes constant as 
our spacecrafts reach farther into interstellar space. Nevertheless, upon a slightly 
more detailed investigation of this behavior, we find that the sinusoidal predictions 
of this model for Pioneer 10 in interval III have amplitude an order of magnitude 
smaller than the quoted value of (0.215 ± 0.022) x 10“^^ km from Ref. 2. Using 
HORIZONS input data for Pioneer 10 from interval III, one finds with a simple non¬ 
linear fit to an annual sinusoid in MATLAB an amplitude prediction of approximately 
0.247 X 10“^^ km s“^. While this most certainly does not rule out the possibility 
that the predictions of this theory can provide a viable alternative to the thermal 
explanation with the ability to model the onset as well as plateau for the anomalous 
‘acceleration’ of Pioneer 11, it does show that we may still need to attribute the 
observed oscillatory behavior to computer mismodeling as was originally suggested 
in Ref. 2. The predicted oscillatory behavior shown in Fig. 2 during interval III for 
which the oscillatory analysis in Ref. 2 was conducted with CHASMP residuals would 
then likely be an order of magnitude too small to discern given the experimental 
precision for Pioneer 10 (see eqn. (52) of Ref. 2). 
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A separate source of concern in this theory’s ability to provide a viable alterna¬ 
tive to the thermal explanation for the Pioneer anomaly comes with its Op predic¬ 
tions for Pioneer 11 located at a heliocentric distance of less than 9 AU. Comparing 
with fig. 7 of Ref. 2, one could initially find discrepancy between the observation of 
almost no anomalous ‘acceleration’ at the first data point near 6 AU versus at least 
50 percent of the total anomalous ‘acceleration’ predicted by the theory of inertial 
centers in Fig. 1(a). However, we must keep in mind the substantial effect that solar 
radiation pressure has on DSN frequency data at 6 AU relative to that of even a 
spacecraft heliocentric distance of 10 AU in our ability to discern a value for Op. 
Referring to fig. 6 of Ref. 2, we realize that, at a distance of 6 AU, the relatively 
large amount of solar radiation pressure could effectively drown out any anomalous 
effect at the level of 10“^^ km s“^ given the experimental precision for Pioneer 11. 
Consequently, the inconsistency between the first data point in fig. 6 of Ref. 2 and 
the theory of inertial centers predicted value in Fig. 1(a) at 6 AU appears to be less 
of a concern than at first glance. 




Fig. 3. Side-by-side plots of the anomalous ‘acceleration’ predicted by the theory of inertial 
centers for Juno (a) and the trajectory of the spacecraft in terms of heliocentric distance (b) from 
6 August 2011 to 1 August 2016. Note that Earth flyby occurred on 9 October 2013.®® 


Focusing our attention on other possible detections of signatures of this theory 
with data from ongoing experiments, we turn to Juno and New Horizons. From 
Fig. 3(a), it is clear that, after Earth flyby on 9 October 2013, this model predicts an 
anomalous jump in ‘acceleration’ for calculated DSN frequency residuals reminiscent 
of what one might expect from an effect resembling the flyby anomaly.^® However, 
before pursuing this idea further, we refer to the amplitude of these changes in 
‘acceleration.’ Given the data intervals used about Earth flyby are on the order of 
days instead of years along with the fact that light-times back and forth to the 
spacecraft are significantly smaller for Juno about Earth flyby than for Pioneer 
10 and Pioneer 11 during the latter portion of their cruise phases out of the Solar 
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System, we find that a change in ‘acceleration’ on the order of 10“^^ km near this 
gravity assist date could not be detected with the fitting procedures employed even 
with a significant improvement in precision from Juno relative to that of Pioneer 10 
and Pioneer 11.®° This reasoning also applies for previous Earth flybys with other 
spacecrafts (see table I in Ref. 35) where similar small changes are predicted by this 
model. It appears at first glance that the theory of inertial centers would be unable 
to explain the flyby anomaly unless some effect in our fitting procedures amplihes 
these small predicted changes in ‘acceleration.’ However, a null signal prediction 
from this theory for Juno Earth flyby would be consistent with the recent analysis 
of Ref. 59 (see Ref. 61 for a separate analysis of the flyby anomaly in the context 
of Juno Earth flyby). 



Fig. 4. Anomalous ‘acceleration’ predicted by the theory of inertial centers for New Horizons 
from 20 January 2006 to 1 September 2015. 


We look to cruise phase after Earth flyby for Juno, where the solar radiation 
pressure begins to decrease in approach to Jupiter. The relevant portion of Fig. 3a 
for us to consider for the predicted anomalous ‘acceleration’ would be from early 
2014 to early 2016. One finds from this figure a potential signature to decipher out 
of the solar radiation pressure between about 2 and 5 AU, referring to Fig. 3(b) 
for Juno heliocentric distance during this time interval. Our ability to discern a 
possible signal of the predictions presented in this article may be suppressed for 
Juno given the magnitude of solar radiation pressure inside of 6 AU. Nevertheless, 
DSN frequency data from cruise phase of New Horizons could provide an additional 
opportunity to test this theory because of the large decrease in solar radiation 
pressure outside of 10 AU relative to the magnitude of the ‘acceleration’ that we 
are searching for as well as the nature of the signal predicted by this model, shown 
in Fig. 4. From a spacecraft heliocentric distance of 10 AU in June 2008 to about 
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33 AU at anticipated Pluto encounter in July 2015, one finds an almost constant 
anomalous ‘acceleration’ prediction nearly of the same magnitude as that observed 
with Pioneer 10. 


Table 2. Replication of table II from the original Pioneer anomaly analysis.^ Error budget for Pioneer 10 and Pioneer 
11. A summary of biases and uncertainties. 


Item 

Description of error budget constituents 

Bias (xlO km s '^) 

Uncertainty (xlO km s 

1 

Systematics generated external to the spacecraft: 




(a) Solar radiation pressure and mass 

+0.03 

±0.01 


(b) Solar wind 


± < 10-5 


(c) Solar corona 


±0.02 


(d) Electro-magnetic Lorentz forces 


± < lo-'^ 


(e) Influence of the Kuiper belt’s gravity 


±0.03 


(f) Influence of the Earth orientation 


±0.001 


(g) Mechanical and phase stability of DSN antennae 


± < 0.001 


(h) Phase stability and clocks 


± < 0.001 


(i) DSN station location 


± < 10-5 


(j) Troposphere and ionosphere 


± < 0.001 

2 

On-board generated systematics: 




(a) Radio beam reaction force 

+ 1.10 

±0.11 


(b) RTG heat reflected off the craft 

-0.55 

±0.55 


(c) Differential emissivity of the RTGs 


±0.85 


(d) Non-isotropic radiative cooling of the spacecraft 


±0.48 


(e) Expelled Helium produced within the RTGs 

+0.15 

±0.16 


(f) Gas leakage 


±0.56 


(g) Variation between spacecraft determinations 

+0.17 

±0.17 

3 

Computational systematics: 




(a) Numerical stability of least-squares estimation 


±0.02 


(b) Accuracy of consistency and model tests 


±0.13 


(c) Mismodeling of maneuvers 


±0.01 


(d) Mismodeling of the solar corona 


±0.02 


(e) Annual/diurnal terms 


±0.32 


Estimate of total bias or error 

+0.90 

±1.33 


3. Conclusions 

To begin our conclusions, we offer remarks on the thermal explanation of the Pi¬ 
oneer anomaly. Refs. 17 and 4 claim that the Pioneer anomaly can be explained 
by a thermal emission term in the line-of-sight acceleration of the spacecrafts that 
decays exponentially and is directed back toward the Solar System. This reasoning 
is consistent with the decay of the onboard ^^®Pu isotope by alpha emission with a 
half-life of 87.7 yr while taking into account the decrease with time in the efficiency 
of converting RTG power to spacecraft power. They conclude the Pioneer anomaly 
is 80 percent thermal based on telemetry data from Pioneer 10 and Pioneer 11 and 
further that there is enough uncertainty in the inference of thermal emission that 
a 100 percent contribution is not unlikely. Furthermore, their logic appears even 
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more convincing after presenting fits to the extended set of Doppler tracking data. 
However, their conclusions are significantly different from those of Refs. 62 and 2 
where thermal emission by the spacecrafts is taken to be 6 percent of the anomaly 
with an uncertainty of 6 percent. The analysis of Ref. 2 is based on blueprints of 
a SNAP-19 RTG power source with symmetrical and larger-than-usual cooling fins 
added specifically for the Pioneer spacecrafts. Each RTG was placed on one of two 
booms extending 3 m from the center of the spacecraft providing a total power 
delivered to the spacecraft bus of 155 W at launch. With a simplified spacecraft 
model. Ref. 2 suggests that a thermal emission of even 60 percent, the one-sigma 
value of Ref. 4, is at the nine-sigma level. Additionally, the simple thermal model 
used in Ref. 2 was applied successfully in quite a different setting to the flyby of 
Gassini about the satellite Rhea,®^ so it does appear likely to be credible. One of 
the authors in Refs. 2 and 62 (ST) worked closely with another author (MMN) 
to develop Table 2, including the thermal model of 6% with a realistic 6% error 
bar. That author later lead the team that developed an alternative model based 
on Pioneer telemetry records. Although we doubt that the telemetry data can be 
stretched beyond a thermal contribution of 70% ± 20%, we nonetheless accept both 
published models as viable. A weighted mean of the two models yields 12.1% ± 
5.7%, which we adopt as the best estimate of the thermal contribution. The effect 
of the much larger estimate of the thermal contribution from the telemetry data 
is to positively bias the earlier estimate in Ref. 2 and Ref. 62, based on a detailed 
model of the Pioneer power subsystem, by one sigma. 

A potential issue with the analysis of Ref. 4 arises if a best-fit inverse square with 
heliocentric distance term is subtracted from their predicted additional acceleration. 
The result of subtracting such a term from their analysis is consistent with the claim 
by Ref. 64 that the anomalous acceleration appears to be nearly constant instead of 
decaying over time. Such a term to subtract away is certainly present from solar 
radiation pressure, but Ref. 4 assume the existing solar pressure during cruise phase 
out of the Solar System for both Pioneer 10 and Pioneer 11 can be taken into account 
by using the model established during the cruise phase of the two spacecrafts between 
Earth and Jupiter. They do not introduce an extra k/d^ term into their fitting model 
for spacecraft heliocentric distances farther than Jupiter, which is a reasonable 
assumption but nonetheless may not be correct. It is important to introduce /c as a 
free parameter and at least consider the possibility that it represents a difference in 
the model for solar radiation pressure between 1 and 5 AU with several parameters 
that could be biased by correlations and an appropriate model for solar radiation 
pressure with only one parameter in the region between 30 and 70 AU. In other 
words, we are not suggesting that Ref. 4 are incorrect in their analysis, but we claim 
instead that their results are not conclusive based on the telemetry data and the 
orbit determinations for the two spacecrafts. Thus, it would seem that the Pioneer 
anomaly remains a mystery until a new deep-space mission makes an unambiguous 
measurement of any extra acceleration or lack thereof, or until a new theory predicts 
the existence of the extra acceleration observed with Pioneer 10 and Pioneer 11 in 
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a manner consistent with all of the observed features of the anomaly. Our purpose 
with this article is to indicate a possible theoretical direction which for now seems 
to be consistent with the observed anomaly. 

We attempted to model the possible onset of the Pioneer anomaly after Pioneer 
11 Saturn encounter using predictions from the theory of inertial centers to offer 
a potential alternative to the thermal explanation of Ref. 4. We derived the form 
of the frequency residuals calculated by the DSN according to this theory and 
discovered a drift term with a clock acceleration similar in magnitude to that of 
Ref. 44 but having an additional weighting factor dependent upon the direction of 
transmission of each photon relative to the direction of our galactic center. However, 
an ambiguity arose with regard to the question of how to relate time as defined by 
the DSN to physical clock-rates in the theory of inertial centers. We adopted a 
particular convention to address the ambiguity that allowed both time coordinates 
to span a “physical range” (i.e. 0 < Tdsn < oo and 0 < t < oo) and, with this 
convention, found that our DSN residual drift term is blueshifted in nature for both 
Pioneer 10 and Pioneer 11 as observed in Ref. 2. After plotting with HORIZONS input 
data, it became clear that our additional weighting factor provided the necessary 
adjustment to produce a jump in anomalous ‘acceleration’ after Pioneer 11 Saturn 
flyby with an eventual plateau causing a constant drift of DSN fractional frequency 
residuals with rate of change similar to the value of the Hubble constant. 

Nevertheless, the additional ‘frequency shift’ according to this cosmological 
model is not in fact a physical shift in frequency of each of these photons. Instead, 
due to the manner in which the DSN calculates changes in frequency, we interpreted 
this predicted shift as a manifestation of the difference in clock-rates between the 
DSN’s definition of time (18) and time according to the theory of inertial centers. 
As a consistency check, we also found our results for Pioneer 10 coincided with the 
observations reported in Ref. 2 as a near-constant clock acceleration for the set of 
data analyzed. Thus, prior to proper fitting of this model to the extended set of 
DSN frequency data mentioned in our introduction, it appears that the predictions 
of the theory of inertial centers could potentially provide an alternative to the ther¬ 
mal explanation of the Pioneer anomaly with the ability to replicate an onset of said 
anomaly after Saturn encounter by Pioneer 11. Still, we are far from satisfied with 
our explanation for why our convention choice somehow resolves the mathematical 
ambiguity one encounters from the restrictions imposed by this theory. In particu¬ 
lar, given that our restrictions on the DSN time coordinate in our latter scenario 
(24) are physically consistent with the expectations from the DSN’s definition of 
range, why should our choice in adopted convention have any physical relevance at 
all? And why can we consider the latter scenario to be unphysical? A completely 
satisfactory resolution of this matter has yet to be determined. 

Future work with regard to the issue of the possible onset of the Pioneer anomaly 
requires fitting of the additional clock acceleration term predicted by the theory of 
inertial centers to the residuals of the extended set of DSN frequency data for 
Pioneer 10 and Pioneer 11. Accurate determination of -x/A should result from this 
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analysis once adjustments to the expected DSN residuals from both spacecrafts are 
incorporated. These adjustments necessary for fitting include the revised thermal 
model for the spacecrafts mentioned earlier in our conclusions with a weighted 
contribution estimate to the anomaly of 12.1% ± 5.7% as well as any adjustments 
to models for the solar radiation pressure such as those suggested in Ref. 64. After 
removing these contributions from the Pioneer anomaly signal in the DSN residuals, 
we should find a fitted value for the Hubble constant according to the theory of 
inertial centers from both Pioneer data sets and subsequently be able to directly 
quantify the goodness of fit of the predictions of the theory of inertial centers to 
actual observations. Additionally, we hope to fit possible releases of DSN frequency 
data from the cruise phases of Juno and New Horizons to the predictions displayed 
in Figs. 3a and 4 to further test this model against observation. 

As remarked upon at the end of Appendix B, a separate line of inquiry into the 
ability of this theory to resolve other observed Solar System anomalies seems to 
appear when incorporating the relative velocity of our observers in DSN-calculated 
fractional frequency residuals. We find another DSN Doppler term that one would 
likely attribute to a fictitious ‘acceleration’ 

^pert — 2|krf|VAvfi, (44) 

where v/{ is the line-of-sight velocity between our observers. This ‘acceleration’ 
term is of a similar form as the empirical solution proposed in Ref. 65 to resolve 
the anomalous secular increases of the astronomical unit and the eccentricity of the 
lunar orbit. These predictions from (44) could become apparent with the analysis 
of Ref. 36 for the secular increase in the astronomical unit as their calculations 
employ significantly more Doppler measurements than the ephemerides, although 
apparently “downweighted to avoid influence of their systematic errors” (see their 
section 3.1). While it is clear from our earlier calculations that the predictions of the 
theory of inertial centers would not affect ranging measurements to the Moon, there 
are subtleties that need to be investigated further with regard to reported increases 
in the eccentricity of the lunar orbit as determinations of lunar GMmoou rely upon 
DSN Doppler data from orbiting spacecrafts.^^’On the other hand, if (44) were 
able to predict the anomalous increases in the AU and lunar orbit eccentricity, why 
then would our larger clock-drift term not be apparent in these analyses? 
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Appendix A. Geodesic solutions 

We look for geodesic solutions in a particular inertial frame of reference as defined 
by the theory of inertial centers.Our geodesic equation is given by 

0 = ^7“VaD^ 


(A.l) 
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where = dx^/da are the components of the affinely parametrized tangent vector 
to our geodesic and x^{a) is our coordinate parametrization of the geodesic. To find 
solutions, we start by expressing this equation in terms of Minkowski coordinates 
(T, X), where we take the spatial coordinate origin of this system to be located at 
the position of our inertial center (i.e. r = 0 —>■ (X, F, Z) = (0, 0,0)). Our metric in 
terms of this Minkowski coordinate parametrization of our inertial reference frame 
in this theory takes the form 

ds^ = -cldT'^ + dX/, (A.2) 


where i = 1,2,3, = X, = Y,X^ = Z, and = —dx^ = —Ar^dr^ for 

massive geodesics. r{T) is the physical distance to our inertial center. As our metric 
components in Minkowski coordinates are all constants, our affine connection tensor 
for Minkowski coordinates is zero: F^^[Minkowski = 0. Our equations of motion in 
Minkowski coordinates are then found to be 



d^X’- 

da^ 


= 0 , 


(A.3) 

(A.4) 


where cr —>■ y = a/A J r(T)dT in the massive case. Examining only spatial compo¬ 
nents, we see that 


which leads us to 


_ d^X^ d UTdX^\ 

da^ da \ da dT ) 

d'^TdX^ fdTVd'^X^ 

^ / dT'^ ’ 

d^X^ 


(A.5) 


(A.6) 


where we have used the equation of motion for our Minkowski time coordinate T. 
The solutions to (A.6) are of the form 


X*(T) = V^T + X/,, 


(A.7) 


where Xg is our spatial Minkowski coordinate location when T = 0 and E® = 
dX"^/dT are independent of T restricted by the normalization condition — 

Cq. Nevertheless, we know from the form of our affine parameter y that we wish 
to express this coordinate parametrization of our geodesics in terms of the more 
physically applicable radial Rindler coordinates where coordinate time t in the ra¬ 
dial Rindler chart (see (11)) progresses at the same rate as a physical clock of 
an observer stationary relative to our inertial center according to this theory (i.e. 
Minkowski coordinate time T does not correspond to a physical clock-rate). The rel¬ 
evant coordinate transformations from our Minkowski coordinate parameterization 












June 29, 2015 0:20 WSPC/INSTRUCTION FILE onsefijmpd 


26 M. R. Feldman & J. D. Anderson 


are 

T = — sinh(-\/At) 
Co 

X = r cosh('\/A<) sin 0 cos 4> 
Y = r cosh(-\/At) sinSsini/) 
Z = r cosh(-\/At) cos 9. 

Using these transformations in (A. 7), we find 

'YA 

Co 


(A.8) 

(A.9) 

(A.IO) 

(A.ll) 


r cosh('\/At) sin0cos0 = ^-J r sinh(-\/At) + Aq (A.12) 

rcosh('\/At) sin0sin()) = fr sinh(-\/At) + Iq (A.13) 

V Co y 

r cosh(-\/At) cos 9 = (| r sinh(-\/At) + Zq. 

V Co y 


(A.14) 


Solving for r by squaring each one of these equations and adding the results together 


' cosh^(-\/At) = ^ ( —rsinh(-\/A<) + A) 


(A.15) 


where we choose our time coordinates T and t to coincide for our initial conditions 
(i.e. t = 0 T = 0 and [dt = dT] |t=r=o)- This requirement forces our Minkowski 
constant co to be equal to where tq is the initial distance of our object to the 

inertial center and is fixed by the value of the Hubble constant. 

We employ the relations 

= Y(^YT)f\T=o = Y(^Yt)n=o = rl (A.16) 

i i i 

__ ^ vi 

^ A*— = ^ AW\ (A.17) 

i i 

where ^^^cl V^iT) = dR/dT. Yet, since dX'^/dT is a constant for all 

T, we have 

^ A* W = RV^\t=o = Rv%=o = rovl (A.18) 

i 

where Vq = dr/dt\t=o is the initial radial velocity of the object in the inertial system 
(i.e. v'‘ = dx^/dt are physical velocities). Plugging in above, we find after solving for 
r{t) 


r{t) = 


Co 


1 + sinh^(yAt)(l - I CO 


— sinh('\/At) 






(«1 

Co 


sinh^(-\/At) + 


1 + sinh' 






(A.19) 
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When t = 0, we want our expression to reduce to tq and not —tq. Therefore, we 
take only the positive root for our solution. 

Additionally, we see that can be expressed in terms of the induced 

Riemannian metric components hij for a space-like hypersurface in our original 
Minkowski coordinates in the following manner: 

= =vl (A.20) 


However, this is simply the squared magnitude of the R®’s in Minkowski coordi¬ 
nates along this hypersurface. Furthermore, since our time coordinates coincide for 
our initial conditions, vq can be interpreted physically as the initial speed of our 
inertially traveling object measured in terms of the clock of a stationary observer 
located at r = tq. In other words, since = dX^/dT = dX^/dT\T=o, 


= E'' 


dX^ dX^ ^ - dX^ dX^ 
"^~dT~dT ~ ^ dT dT 


= EE^ 


dX^ dx^ dX^ dx^ 


i,j k,l 

= EfE^ 


dx>^ dT dx^ dT 


kA 




dX^ dX^ \ dx^ dx‘ 

)~dr~di 


T=0 


T=0 




E dx^ dx^ ^T 

'‘’■•wit , =L''«» 

k,l *=0 




= (^'o)^ + ^0 + sin^ , (A.21) 

where hyi are induced Riemannian metric components along the t = T = 0 hy¬ 
persurface for the radial Rindler chart and x^ are spatial coordinates in the same 
chart. Plugging back into our expression for r{t) where we remember to take only 
the positive root, we find for the radial position of our object 

ro 


r{t) = 


1 -I- sinh^(-\/At)(l — (^)^) I Co 


• Sin 


h(-\/At) 


M + sinh^(-\/At) 


I 


(A.22) 


. Co y V Co, 

with Co = VAxo- 

As briefly alluded to earlier in this section, our normalization condition for the 
tangent vector to a geodesic is 

-k = gabU‘^U\ (A.23) 

where 

0 null geodesics 
1 time-like geodesics 


k = 


(A.24) 
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We can express this normalization condition in our original Minkowski coordinates 
in the following manner: 

2 / \ 2 


1. 2/^ 


E 


/dT 
\ da 




EV 

da ) 


dX^ 

da 

Co - Vo 


which provides us with the constraints 

Vo = Co massless objects 
Vo < Co massive objects. 


(A.25) 


(A.26) 

(A.27) 


Returning to our coordinate transformations, the angular positions of our object as 
a function of inertial time t are found to be 


(j){t) = tan 


^r(t) sinh(yAt) + Yo 


= 00 + tan 


-1 


-^r(t) sinh(-\/At) + A^o- 

<p 

r sinh(vAt)^co sin 6*0 

ro sin 00 + r sinh(-\/At)[^ sin0o + cos0o] ^ 


and 


0{t) = cos 


-1 


■^0 


— tanh(-\/At) H- ^ 

Co r(t) cosh(vAt)_ 


= cos 


-1 


ro cos 00 


/u’' v^ \ 

tanh(-\/At)( — cos0o-2.7.psj]2 0 Q ) + , 

yco Co / rcosh(vAt)J 


where we’ve employed 


X (t) |t=o = Ao = ro sin 0o cos 0o 

Y{t)\t=o =Yo = rosin 00 sin 00 
Z{t)\t=o = Zo = ro cos 00 , 


and 


= Vq = Vq sin 00 cos 0 o + r^ro cos 0 o cos 0 o — VqVo sin 0 o sin 0 o 


V = Vq = Vq sin 00 sin 0o + VqTo cos Oq sin 0o + VqTo sin 0o cos 0o 

= Vq = Vq cos 00 — rQrosin0o. 


(A.28) 


(A.29) 

(A.30) 

(A.31) 

(A.32) 

(A.33) 

(A.34) 

(A.35) 


Appendix B. DSN frequency shift with the Doppler effect 

To factor in the Doppler shift from the relative velocity of our observers, we return 
to the general expression for the observed wavelength of a photon according to the 
theory of inertial centers taken in Ref. 43 to be 

27r 
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where C/“ is the ‘four-velocity’ of the observer affinely paramterized by y (not simply 
t) and given by = dx^/dx in component form, is the wave-vector of the 
photon, and P refers to the point of emission/reception of the electromagnetic 
signal in question. We then have the following relations as the observer is assumed 
to follow a time-like geodesic while the photon propagates along a null geodesic: 


C/“17a = -1 
k^ka = 0 . 


Expanding these relations in the radial Rindler chart, one finds 


W = 




Kr^ — cosh^('\/At) ( + sin^ 9{v'^) 


1 - 1/2 


(B.2) 

(B.3) 


(B.4) 


Ar^ = {Pf + cosh2(VAt) -h sin^ , (B.5) 

where we define u* = dx''/dt and c® = employing the Latin index i in the 

component form of these vectors to refer to spatial components. Notice that is 
the local velocity of our observer relative to the inertial center about which he/she 
moves (the center of the Milky Way for our purposes). Yet, along a constant t spatial 
hypersurface with induced Riemannian metric components hij , one finds 

* 

= (c’’)^ -I- cosh^(-\/At) ^(c®)^ -I- sin^ , (B.6) 

and thus 


J2c"c^ = Ar^. (B.7) 

i 

Note that c here is not the DSN-adopted value for the speed of light in empty 
vacuum. Instead c in this section refers to the speed of the photon at a particular 
space-time point in the inertial frame according to this theory. For our own ease, 
we also denote the square of the local velocity of our observer relative to the center 
of the inertial system by 


^ v^Vi = + r^ cosh^('\/At) f (u®)^ -I- sin^ 9{v^Y 

i ^ 

Then our expression for t/* can be written in a simpler form. 


[/* = 


^fkr 

and the wavelength measured by our observer is found to be 


2tt 


•\/Ar 


A Y^I — 


1 - 




(B.8) 


(B.9) 


(B.IO) 
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where = Si j hijv^c^ can be thought of as the ‘dot-product’ of the observer 

and photon velocities relative to the center of the inertial system along the constant 
t spatial hypersurface. To find an expression for k*, we make use of the Killing 
vector field = (l/'/Xr) cosh{'/At){d/dt)°'—smh{'/At){d/dr)°' in the conservation 
equation —E = ^‘^ka, where E is a constant."'’^ One finds 




E 

\/Ar[cosh(-\/At) -I- (c^/VAr) sinh(-\/At)] ’ 


(B.ll) 


and our relation for the observed wavelength reduces to 


27r 

T 


E 


cosh(-\/At) J- (c^/VAr) sinh(-\/At) ^1 — v^/c^ 


1 - 


Si ^'*Ci 


(B.12) 


We recall U* = dt/dx = (I/a/Ac) • dt/dr, where r is the proper time according 
to the theory of inertial centers. Comparing this expression for U* with (B.9) and 
plugging into (B.12), we have 


27r 

T 


E 


dt 


1 - 


S, 


U‘Ci 


cosh(-\/At) + {c^/\/Ar) sinh(-\/At) 

Taking the ratio of observable wavelengths, one finds for the uplink shift 

cosh(-\/Ati) + {cl/\/Ari) sinh(-\/Ati) {dt/dT )2 


Ai 

A 2 


uplink cosh(A/At 2 ) + (c2/\/Ar2)sinh(A/At2) (dt/dT)i 


[1 



-1 

L C2 J 

2 

L C2 J 



and for the downlink 
A 2 

At 

downlink 


cosh(-\/At 2 ) — (c 2 /\/Ar 2 ) sinh(-\/At 2 ) {dt/dr)^ 


cosh(-\/At 3 ) -b (cg/v^ca) sinh(-\/At 3 ) {dt/dT )2 


1 - 


s. 


3 L 


1 -b 


s. 


(B.13) 


(B.14) 


(B.15) 


where we have implicitly taken c^| downlink = -C2I uplink and I uplink = I downlink- 
Furthermore, all quantities marked with a ‘2’ in the rest of our analysis refer to 
the uplink portion of the two-way traversal, and we assume that ki = k 2 |upiink = 
—k 2 1 downlink = —kg, u*|i = u^jg, and rg = rg which implies ci = 03 . Our expression 
for the two-way shift then reduces to 


Ai 

A3 


cosh(-\/Ati) -b (Cd/-\/Ari) sinh(-\/Ati) 
cosh(-\/At 3 ) — (ci/\/Ari) sinh(-\/At 3 ) 

(dt/dr)3 

{dt/dr) I 


cosh(-\/At 2 ) — {c 2 /y/Ar 2 )sYaE{y/At 2 ) 
cosh(-\/At 2 ) + (c 2 /\/Ar 2 ) sinh(-\/At 2 ) 


1 + (Si^'^c^/c^)! 1 - (S* 


1 - (Si v'-Ci/c^)i 1 -b (S* vVi/b^2 


^ 6 ) 
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Simplifying the first line of the above expression with Taylor series expansions up 
to first-order in '/Kt, we have 


^3 


1 -I- (kl • f)V^[(t3 - t2) - {t2 - tl)] + 0[M^) 


{dt/dT)3 

{dt/dT)i 


1 - 1 + (Ei ^^*C*/c2)2 ’ 


(B.17) 


where we have substituted for = Ci(ki-f) = •\/A?’i(ki-f) and C 2 = 

C2(k2|upiink-f) = A/Ai’2(k2 1 uplink-f) = \/Ar2 (ki-f) , taking the galactic radial unit 
vector to be approximately the same for our observers in the Solar System (i.e. 
fi = f 2 = f)- Yet from our earlier analysis in the discussion of this article, we 
know that the elapsed coordinate time for the uplink, tupiink, will be approximately 
equal to the elapsed coordinate time for the downlink, tdowniink- In other words, 
ts —t 2 = t 2 — tl resulting in 


Ai idt/dT)3 1 + 1 - (E»'^'c,/c^)2 .O.A.2N 

A3 {dt/dT)i ■ 1 - (E, v^c,/c^)i ■ 1 + (E. v^cjc^)2 ^ ^ 


(B.18) 


However, for relative frequency shifts calculated by the DSN, one must realize 
that the DSN defines physical clock-rates from relation (18). Thus, while our actual 
fractional change in frequency once formulating in terms of coordinate time cor¬ 
responding with physical stationary clock-rates according to the theory of inertial 
centers will not produce a drift term as shown later in this section, what the DSN 
calculates in terms of the clock-rates it assumes to be physical will. In other words, 
the range-rate calculated by the DSN uses the DSN’s definition of time for its time 
coordinate. However, the DSN’s time coordinate drifts away from actual physical 
time according to the theory of inertial centers, and thus with this DSN calculation 
we no longer have an accurate assessment of the frequency observable itself (i.e. 
vzjvi 7 ^ [r's/j^iJIoSN)- To understand what the DSN calculates, it appears we must 
substitute t —>■ Tdsn in our previous expression as T = Tdsn defined by (18) is 
the coordinate time that the DSN uses in its range-rate calculation. Still, in our 
previous expression, r remains the physical proper time according to the theory of 
inertial centers. We also keep in mind that u* = dx'^/dt is the local velocity of the 
observer in the inertial frame. This local velocity happens to coincide with the ve¬ 
locity expected from the DSN when examining for our slow-moving massive objects 
(i.e. dr/dr Ri 0 in this limit in (9)). Then, 


Ai 


A3 


DSN 


{dT/dT)3 ^ 1 + iT,^'^’'ci/c^)l _ 1 - (Ew'd/^^)2 
{dT/dT)i 1 - (Ei'^*Ci/c2)i 1 + (Ei^*D/c2)2 


C>(At2).(B.19) 


If we keep terms up to first-order in (f/c), we see from our earlier expression for {7* 
in (B.9) along with the form of our affine parameter in ( 6 ) that 


dt 

dr 


i+o[{v/cn 


(B. 20 ) 
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Thus, 
Ai 

As 


{dT/dt)3 I + 1 - 


0[Kt\{vlcl%2l) 


iJDSN {dT/dt)i I - l + (Ei^^*Ci/c2)2 

Finally, we use our expression for the difference in clock-rates between time 
according to the theory of inertial centers and time assumed by the DSN. From 
(30), we have 


dT'\ 
dt J, 


VArs 

CDSN 


1 -I- \/A|ki-r|tj 


OiAf), 


(B.22) 


where we adopt our former convention of = constant. Furthermore, for the uplink. 


f ri r 2 > ri 

\ r 2 ^2 < ri ’ 


(B.23) 


and the analogous situation with results for the downlink, although one need not 
worry too much about this since these constant terms cancel in our two-way DSN 
frequency shift expression. Plugging in above and Taylor expanding our fractions 
to first-order in (u/c) and VAt, 


Ai 


.As. 

DSN 


1 -I- A/A|krf|(t3 - ti) 


1 + 2(^?;*Ci/c^)|i - 
i 


1 -I- A/A|krf|(t3 - ti) 


1 -b 2(vi-ki/ci) - 2 (v 2 -ki/c 2 ) 


2(^u*c*/c2)|2 

i -* 

+0[Af, (u/c)2] 

+ 0[AtMv/cm.2A) 


where we have taken the spatial hypersurface of the metric in our local region of 
the Milky Way within our solar system to be approximately Minkowskian resulting 
in being equal to the ordinary Euclidean dot product v-ck. In other words, 

the difference in expectation for what constitutes a physical clock-rate does not 
affect the Doppler terms to first-order in '/At. To see why we can approximate 
in this manner, we notice that the difference between the spatial portion of the 
Minkowski metric transformed into spherical coordinates and that of (II) lies in 
the cosh(-\/At) term. Yet, since we are examining up to only first-order in /At 
after replacing t ^ T for DSN frequency calculations, these two different spatial 
hypersurfaces will appear to be the same. Furthermore, since all of our calculations 
are for spacecrafts traveling near the Sun relative to the distance to our galactic 
center, we take (v/ci) ~ {v/c/) as we’ll have ri/r 2 multiplying factors of v/c in this 
expression to relate terms with ci = /Ari to those with C 2 = VAr 2 . Since (v/c) 
is already small for our slow-moving massive objects, this approximation appears 
valid. Our general expression for the DSN-calculated two-way shift in wavelength 
according to the theory of inertial centers is given by 


= 1 - 


2 kr(v 2 - vi) 


DSN *^1 

4A/A|ki-r|kr(v2 - Vi) 


(O - 


f 2V^|ki-f|(t2 - ti) 
h)+0[AtMv/c)% 


Cl 


(B.25) 
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where we have used the approximation — ti = 2(^2 — ti) for electromagnetic sig¬ 
nals traveling within our localized region of the Milky Way. For this DSN two-way 
wavelength shift, we find our stationary clock-drift which we derived in our discus¬ 
sion, 2-\/A|ki-f I (^2 — t i), in addition to the DSN-assumed expression [Ai/AsJdsn = 
1 — 2ki-(v2 — vi)/ci, as well as another term —4-\/A|ki-f|ki-(v 2 — Vi)(t 2 — ti)/ci 
dependent upon the product of the line-of-sight velocity of the spacecraft and our 
stationary clock-drift term. 

It is interesting to note that if we were to interpret the cause of this last term 
in the DSN Doppler shift as the result of an additional relative ‘acceleration’ acting 
between our two observers, we would likely assume this ‘acceleration’ to take the 
form 

^pert — 2 |krf|VAvfl, (B.26) 


where vr is the line-of-sight velocity between our observers. Here, we’ve used the 
assumed form of residual DSN Doppler terms due to an unmodeled acceleration, 

^pert; 

= (B.27) 

model ^ 

where R = ki is the line-of-sight unit vector. Notice that (B.26) is of a similar 
form as the empirical solution proposed in Ref. 65 to resolve the anomalous secular 
increases of the astronomical unit and the eccentricity of the lunar orbit 


observed 


^lorio — 


(B.28) 


where, in our case, k = 2|ki-f| is similar in value to lorio’s empirical solution of 
2.5 < fc < 5. Of course, in our notation, '/A = Hq where Hq is the Hubble constant. 

Furthermore, we wish to emphasize that this ‘frequency shift’ is the result of the 
DSN coordinate parametrization of time and not a physical shift in frequency. To 
see this, we remark that the calculation for the actual shift in frequency according 
to this theory requires that we use coordinate time t if we assume that coordinate 
time corresponds with the physical clock of a stationary observer. Returning to our 
general expression (B.18) for the fractional shift in wavelength according to the 
theory of inertial centers, we realize that coordinate time t for the initial emitter 
moving with the same velocity v'^\i = and remaining the same distance away 
from an inertial center ri = has the same relation with proper time r for emission 
and reception such that (dt/dr)|i = {dt/dT)\^. Then, to first-order in and (u/c), 
one has for the two-way fractional shift in wavelength after factoring in the relative 
velocities of our two observers. 





+ 0[At^{v/c)% 


(B.29) 


where we obtain the DSN-assumed expression for ri k, r 2 . 
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